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Using  an  analytical  model  to  represent  the  soil  and  an  clastic- 
membrane  theory  to  represent  the  sheath  commonly  used  around  the  soil 
sample  the  stability  of  the  soil-membrane  system  is  studied.  The  measure 
of  stability  used  is  the  time  derivative  of  the  rate  of  doing  work  by  the 
external  loading  for  a  prescribed  motion. 

The  results  give  a  means  of  evaluating  the  relative  effect  of  the 
membrane  on  the  initial  stability  of  the  ststem.  A  numerical  example  shows 
the  effect  tc  be  negligible  in  the  case  of  the  initial  stability  of  a  par¬ 
ticular  soil, 
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Introduction 

In  an  earlier  work  [i]*  the  authors  have  attempted  to  generalize  the 
usual  approach  to  time  independent  continuum  soil  mechanics.  It  has  been 
popular  recently,  for  example  see  [2],  to  use  a  perfectly  plastic  Coulomb- 
Mohr  theory  with  an  associated  flow  rule  determined  from  the  normality  of 
the  incremental  plastic  strain  vector  to  the  yield  surface.  The  theory 
proposed  by  the  authors  differs  from  such  an  approach  in  the  following 
respects 

a.  The  loading  surface,  instead  of  diverging  with  increasing 
hydrostatic  pressure,  closes  and  is  intersected  by  the 
spherical  stress  axis  at  specific  value  of  hydrostatic 
pressure.  Thus  the  loading  surface  is  a  closed  surface. 

b.  The  size  of  the  loading  surface  varies  for  each  point 
according  to  the  value  of  the  specific  volume  at  that 
point . 

c.  The  normality  of  the  incremental  strain  vector  to  the 
loading  surface  at  each  point  is  replaced  by  a  set  of 
permissible  directions  which  happen  to  always  include 
normality.  The  permissible  directions  are  deduced  from 
energy  considerations  that  include  an  internal  energy  U 
which  depends  on  the  specific  volume. 


*  Numbers  in  brackets  designate  references  listed  in  the  bibliography. 
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U'ing  a  specific  case  of  the  above  theory  the  deformations  occuring  in 

a  triaxial  test  (a  compression  test  in  the  presence  of  hydrostatic  pressure) 

were  studied  [3]  by  the  authors.  This  work  neglected  the  presence  of  a 

membrane  around  the  oylindrical  sample.  In  this  study  of  the  triaxial  test 

two  possible  deformation  modes  were  considered.  The  first  was  a  homogeneous 

deformation  of  the  entire  sample ,  shown  in  Figure  la,  while  the  second  was 

a  shearing  motion  of  an  inclined  layer  of  the  material,  depicted  in  Figure  lb. 

The  results  of  sample  calculations  were  given  for  an  analytical  model  whose 

properties  were  similar  to  a  uniform  Leighton  Buzzard  sand  which  has  been  re* 

ported  in  [4].  For  the  model  the  homogeneous  deformation  indicated  an  initially 

positive  value  for  the  time  rate  of  change  of  the  rats  at  which  the  external 

loads  do  work  on  the  soil  sample.  For  the  same  quantity,  which  will  be  denoted 
•• 

by  W,  the  shearing  motion  indicated  an  initially  negative  value.  A  positive 
•  • 

value  of  W  indicates  a  stable  system  while  an  unstable  system  is  predicted  by 

•*  »« 

a  negative  value  of  W.  The  actual  magnitudes  of  W  serve  also  as  measures  of 
the  stability  or  instability  of  the  system.  Consequently  of  all  possible 
modes  of  motion,  the  mode  with  the  lowest  measure  of  stability,  H,  is  hypo¬ 
thesized  as  the  one  which  in  actuality  will  occur.  For  the  shearing  motion, W 
was  also  shown  in  [3]  to  be  inversely  proportional  to  the  thickness  H  of  the 
shearing  layer  as  shown  in  Figure  lb.  For  the  numerical  example  where  W  was 
shown  to  have  a  negative  sign,  the  value  of  H  producing  the  largest  negative 
W  is  H  equal  to  zero  leading  to  W  equal  to  minus  infinity.  Thexefore  for  the 
modes  considered  a  displacement  discontinuity  across  an  inclined  plane,  as 
depicted  in  Figure  2,  is  the  anticipated  actual  mode  of  deformation. 

In  connection  with  the  shearing  mode  of  deformation  a  change  in  loading 
procedure  develops  which  is  not  present  ir  th<  homogeneous  deformation.  The 
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The  initial  state  of  loading  is  a  vniform  hydrostatic,  -o^,  and  a  superimposed 
compressive  stress,  -c^,,  for  all  modes  of  deformation.  In  order  to  induce 
the  shearing  mode  a  non-trivial  rate  of  change  of  shearing  stress  °yZ*  see 
Figure  2,  must  be  present.  It  is  assumed  in  the  following  that  this  loading 
is  supplied  by  the  testing  device  while  the  sample  undergoes  the  initiation 
of  shearing  motion. 

It  is  the  purpose  of  the  current  work  to  present  an  evaluation  of  the 
effect  of  a  membrane  surrounding  ■"he  sample  in  the  case  of  the  shearing  motion 
in  the  triaxial  test.  The  impermeable  membrane  is  generally  employed  to  help 
maintain  the  integrity  of  the  soil  sample  before  loading  commences  and  to 
sustain  the  hydrostatic  loading  applied  to  the  sample.  That  is,  in  the  tri¬ 
axial  test  there  is  a  pore  pressure  in  the  sample  which  is  maintained  under 
prescribed  conditions  (in  this  work  the  pore  pressure  is  always  zero)  while  the 
external  hydrostatic  pressure  is  applied  to  the  sample  independent  of  the  pore 
pressure.  In  the  current  considerations  both  the  pjre  pressure  and  the  external 
hydrostatic  pressures  are  held  constant  so  that 


hydrostatic  pressure  -  pore  pressure  =  -Oj 
throughout  the  test. 


(1) 


Time  Independent  Soil  Under  Consideration 

The  description  of  the  soil  under  consideration  is  reviewed  in  this 
section.  More  detail  concerning  this  description  may  be  found  in  [3].  The 
functions  used  are  as  follows: 

v  =  specific  volume,  volume  per  unit  weight  of  soil  sample. 

♦  s  4(v)  s  friction  angle  for  Coulomb-Mohr  part  of  loading  surface. 


c  »  c(v)  *  cohesive  strength  for  Coulomb-Mohr  part  of  yield  surface. 

oq  *  ®0(v)  *  mean  normal  stress  associated  with  the  closing  of  the 

yield  surface  on  the  compression  side. 

U  a  u(v)  a  internal  energy  per  unit  weight  of  the  soil  sample. 

o  a  1/3  (o  ♦  o.  ♦  o  )  a  mean  normal  stress 
>?  y  z 

6  «  G(v,o)  a  function  which  determines  the  direction  for  the  incremental 

strain.  When  G  is  zero  the  vector  is  normal  to  the  loading  suriace. 

When  G  equals  2  sin}/[3  (1  +  sin#)]  the  flow  is  incompressible. 

From  [3],  0  i  G  a  a  (3U/3v)  (2  sin#)  /  (3  {1  ♦  sin#}*{o  -  o  t  a 

o 

[3U/3v]})  $  (2  sin#)/(3[l  +  sin#])  and  a  is  a  constant  for  this  case. 
Z,  R,  H  a  dimensions  shown  in  Figure  lb. 

Z  a  constant  vertical  velocity  of  compression  loading  device ,  Figure  2. 

Y  a  lateral  velocity  of  section  of  sample  above  shearing  plane, 

Figure  2. 

The  following  conditions  are  assumed  satisfied  in  this  investigation. 


c(v)  i  0 

(2. a) 

0  S  #  (v)  <  w/2 

(2.b) 

3v  *  ° 

(2.c) 

3U 

O  <  -r-  <  0 

o  3v 

( 2 .  d  ) 

o  >  a 

o 

(2.e) 

0  <  a  <  1 

(2.f ) 

1  >  ic  >  0 

(2.g) 

In  [3]  it  is  shown  that 


where 


S**  =  (Z/H)  [ oj  (Z/H)  +  {o  /(Z/H)>  {(Y/H>/(Z/H)}] 
L  yz 


(4) 


and 


°c  *  -2lc  -  °L> 


o  =  oc/3  +  oL 


sin$ 


Jl  -  sin$ 


*  ■  -  -1  ou!  ;in:r-> 


v/(z/H)  = 

COS^l 


(Y/H)/(Z/H)  -  -tan  * 

Ln$ 
■sin$ 


\/dm  ■•atpSfsrfc +  «-«i>  7~f7J^ 

(l-sin$) 


(5. a) 

(5.b) 

(5.c) 

(5.d) 
(5.e) 
( 5.f ) 


o  /  (Z/H)  r  (0_/2  sin24») 

yz  l 


[cos* 


,OG  •  2  t  G  cos  2»)  H  •  • 


2  -  G  (l  +  sin$)  3v 


(v/CZ/H]) 


.  (1  +  sin$)(3  +  cos2*>  f3G  .  1  3G  ,•  ,r.i  ,r 

*  ~"-sTi-rs'iiif;  (s7 <v/|:i/H])  +  3U(,,c/tz/rtl)l1  (5-g> 

Consequently  when  the  functions  $(v),  c(v),  oq(v),  U(v),  and  G(v,o)  are 
known  for  the  material  then  specification  of  v  and  combined  with  the  above 
equations  determines  (H/Z)  S**. 

In  [3]  the  following  functions  were  used  for  the  sample  calculation 
where  v  is  in  cubic  inches  per  pound: 

$(v)  -  0.686  -  0.100  (v  -  16.08)  {radians} 
c(v)  =  0 


°o(v)  ■  -129(^l6?0?)  tlbs/1"2) 


(6. a) 
(6.b) 
(6.c) 


oU(v)  =  +12.9[-(19.00  -  v)  +  2.92  ln{2.92/(v-16.08) }]  (inches)  (6.d) 
3U  . _ 2sin4» _ 


G(v,o)  * 


3(1  +  sin$)(o  -  o  +  a-|^0 

O  iV 


(6.e) 
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with  the  following  conditions 
v  *  16.83  in3/lb 

°L  *  •50  lb>/in*  <?.b> 
These  functions  and  conditions  led  to  the  following  results  for  the  shearing 
moda  of  deformation 


(£)  S**  =  -1023 

(8. a) 

+  *  34.7° 

(8,b.' 

0Q  *  -343  lbs/in2 

(8.c) 

oc  =  -1J2  lbs/in2 

(8.d) 

o  *  -94  lbs/in2 

(8.e) 

*  =  60.4° 

(8.f) 

v/(Z/H)  *  -17.48  in3/lb 

(  8 .  g ) 

(Y/H)/(Z/H)  s  -1.755 

(9.h) 

aQ/  (Z/H)  =  -774  lbs/in2 

(8.i) 

o  /(Z/H)  *  +  142  lbs/in2 
y* 

(8.  j  ) 

Inclusion  of  Membrane  Forces 

In  this  analysis  the  membrane  is  treated  as  a  thin  walled  elastic 
cylinder.  For  simplicity  in  describing  the  membrane  behavior  consider  the 
single  plane  element  shown  in  Figure  3.  If  the  bending  stiffness  of  this 
element  ic  neglected  then  the  deformation  elements  6£,  6^,  and  6b  shown  in 
Figure  3a  induce  the  following  forces  shown  in  Figure  3b. 


N  *  6£  S  H  E/H 
T  *  6^  S  h  G/H 


where  h  is  the  membrane  thickness,  E  is  the  membrane  Young's  modulus  and 
G  is  the  membrane  shear  modulus. 


(9. a) 
(9.b) 
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The  motion  depicted  in  Figure  lb  may  be  expressed  in  terms  of  Y  and 
Z.  The  displacement  rate  of  the  upper  section  of  the  cylinder  denoted  by 
t  is 

t  *  (o)  i  +  (y)  J  +  (z)  t  do) 

It  is  useful  to  define  a  set  of  orothogonal  unit  vectors  n,  t,  and 
e  as  shown  in  Figure  4.  Set 

n  *  (cos0)  t  +  (sin6)^  +  (0)  k  (11. a) 

t  *  [cos2*  *  sin^cos2©]”1^2  [(-sinOcos*)  1 

+  ( cos 0 cos*)  J  +  (-cosOsin*)  )c  ]  (11. b) 

4  o  2  2  -1/2  -*■ 

e  =  [cos  *  +  sin  *  cos  0]  [(-sin0cos0sin*)  I 

+  (cos  0sin*)  j  +  (cos*)  k  3  (11. c) 

At  the  instant  flow  begins  the  deformation  rate  components  6  and 
„  • 

6^  may  be  defined  as  the  projections  of  ?  onto  the  e  and  t  unit  vectors 
respectively.  Hence 

*  £cos2^  +  3in2^cos20]”^^2  »[Ycos20sin*  +  Jcos*)  (12. a) 

s  s  [cos2*  +  sin2*cos203  ^2  *[Yc(.s0cos*  -  Zcos0sin*)  (12. b) 

Introducing  a  unit  vector  normal  to  the  plane  of  shear 
p  *  (0)  i  +  (sin*)  *  t  (cos*)  ic  (13)  . 

a  relationship  between  H  in  Equations  (9)  and  the  thickness  of  the  shearing 
layer,  H ,  may  be  determined 

H  =  H  e»p  *  H  (sin2*cos*0  +  cos2*)^2 


(14) 
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Using  Equations 


6 

£ 


(9)  and  (14),  and  5  may  also  be  defined  as 


*(cos20sin2*  +  cosV1/2 
hE 

'\cos20sin2*  +  cos2*)”^2 
hG 


(15. a) 

(15. b ) 


Tha  increment  of  length  a\ong  the  ellipse  formed  by  the  intersection  of  the 
shearing  plane  and  the  cylinder  may  be  expressed  in  terms  of  e  by 

dS  s  R(cos28sln2*  »  cos2*)*^2 

cos*  *de  (16) 


The  rate  of  change  of  the  force  required  to  deform  the  membrane  is  denoted 


by  P  and  is  defined  by 


F  = 


9=0  *  +  (^>  ?  J  dS 


(17) 


Making  use  of  Equations  (11),  (12),  (15)  and  (16),  Equation  (17)  may  be 
evaluated  to  yield 


F  =  ifc os?  C{0)*  +  {*3/4)  E  sin2'!'  +  G  cos2*)  Y  +  (E  -  G)(cos*sin*)Z)  } 


+  ((E  -  G)(cos*sir.*)Y  +  (2Ecos2*  +  Gsin2*)Z)k] 


(18) 


Using  Equations  (10)  and  (18),  the  time  rate  of  change  of  the  rate  at  which 
the  external  loads  do  work  on  the  membrane  may  be  evaluted.  Denoting  this 
quantity  as  W  membrane  it  can  be  shown  that 

W  s  p.j 

membrane  ° 

a  *R2Z(Z/Z)2(Z/H).(h/R)  [(Y/z)2  3Esin2»  +  4 

4  COS*  + 

2  (Y/Z)  (E  -  G)  sin*  +  ffi508**  +  Osin2*  , 

cos*  J 


(19) 
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The  quantity  treasuring  the  stability  of  the  soil-membrane  system  is  simply 

th.  su*  of  w50il  and  V  MBhran,.  Thu. 

•*  ••  •• 

U  S  u  +  u 

"soil  membrane 

=  «RJZ  <Z/Z>2  [S«  +  (Z/H)(h/.R)  (CY/Z)2  •  ^*t°  j^*"*”*  *  ♦ 

*  « -  sw 

Using  the  numerical  results  of  the  example  cited  in  the  previous  section 
with 

E  *  200  psi  (21. a) 

G  =  66.7  psi  (21. b) 

yields 

.  n  .  o  u 

(22) 


W  *  vR2Z(Z/Z)2  (Z/H)  (-1023  +  702  |>) 


For  cases  of  practical  significance  a  typical  value  of  h/R  is  0.014. 
Clearly  any  h/R  of  interest  is  too  small  to  make  the  sign  of  (-1023  +  702  h/R) 
positive  and  thus  stabilize  the  system.  The  result  is  that  any  typical 
membrane  has  a  very  small  stabilizing  effect  on  the  initial  stability  of  a 
triaxial  test  for  the  material  and  initial  state  considered. 
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FIGURE  CAPTIONS 


Figure  1  -  Two  Modes  of  Deformation  for  the  Triaxial  Test: 

a)  Homogeneous  Deformation 

b)  Shearing  Motion 

Figure  2  -  Mode  of  Deformation  for  Shearing  Motion  When  the  Corresponding 

«• 

W  is  less  than  zero. 

Figure  3  -  Plane  Element  of  Membrane  Showing 

a)  Displacements 

b)  Forces 

Figure  4  -  Section  of  Thin  Membrane  Below  Shearing  Plane  Showing  Vectors 


BEFORE 


f 


a)  HOMOGENEOUS  DEFORMATION 


BEFORE  AFTER 

b.)  SHEARING  DEFORMATION 


FIG.  I  TWO  MODES  OF  DEFORMATION  FOR  THE 
TRIAXIAL  TEST. 


FIG.  2  MODE  OF  DEFORMATION  FOR  SHEARING 

MOTION  WHEN  THE  CORRESPONDING  W<0 


h—  S 


o.  b. 

FIG.  3  PLANE  ELEMENT  OF  MEMBRANE  SHOWING 

a. )  DISPLACEMENTS 

b. )  FORCES 


n  NORMAL  TO  CYLINDER 
SURFACE 

t  //  PLANE  OF  SHEARING 
AND  «L  rf 

n,7,'e  ARE  EACH  UNIT  LENGTH 

AND  MUTUALLY  PERPENDICULAR 


FIG. 4  SECTION  OF  THIN  MEMBRANE  BELOW 
^SHEARING^PLANE  SHOWING  VECTORS 

n,  t,  AND 


